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Abstract. —The white Barn Owl subspecies (Tyto alba alba) is found in southern Europe and the reddish- 
brown subspecies ( T. a. guttata ) in northern and eastern Europe. In central Europe, the two subspecies 
interbreed producing a large range of phenotypic variants. Because of the different ratios of the sub¬ 
species in different geographic regions, we predict that genetic variation should be greater in Switzer¬ 
land than in Hungary. We tested this hypothesis by measuring genetic variation with the RAPD method. 

As predicted, the genetic differentiation within a Swiss population of Barn Owls was significantly greater 
than the variation within a Hungarian population. This suggests that gene flow is greater in central 
Europe than at the eastern limit of the Barn Owl distribution in Hungary. In both countries genetic 
variation was more pronounced in females than in males. As in other birds, this is probably because 
female Barn Owls are less philopatric than males. The number of migrants between Hungary and 
Switzerland is ca. 1 individual per generation; if calculated separately for the sexes, then 0.525 for males 
and ca. 1 for females (Nm values). The difference in the number of migrants between genders again 
is likely a consequence of higher male philopatry. The sexual differentiation is greater in the Swiss 
population than in the Hungarian and the genetic substructuring of the populations of the species is 
substantial. The reason for the considerable population substructuring could be the nonmigratory be¬ 
havior and socially monogamous pairing of the species, as well as the geographical barriers (Alps) 
between the populations examined. 
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VARIABILIDAD GENETICA ESTIMADA MEDIANTE RAPD ENTRE SEXOS Y POBLACIONES DE 
TYTO ALBA EN EUROPA 

RESUMEN.—La subespecie Tyto alba alba se encuentra en el sur de Europa, y T. a. guttata en el norte y 
oriente de este continente. En el centro de Europa, las dos subespecies se entrecruzan y producen una 
amplia gama de variantes fenotipicas. Debido a la variacion geografica en la proportion de individuos 
pertenecientes a las distintas subespecies, predijimos que la variabilidad genetica deberia ser mayor en 
Suiza que en Hungria. Pusimos a prueba esta hipotesis midiendo la diferenciacion genetica usando el 


1 Email address: bobmatix@freemail.hu 


142 




June 2005 


Barn Owl, Genetic Variability 


143 


metodo de ADN polimorfico amplificado aleatoriamente (RAPD, por sus siglas en ingles). Como lo 
habfamos predicho, la diferenciacion genetica dentro de una poblacion suiza de T. alba fue significati- 
vamente mayor que la variacion observada dentro de una poblacion hungara. Esto sugiere que el flujo 
genetico es mayor en el centro de Europa que en el limite oriental de la distribution de T. alba en 
Hungria. En ambos pafses, la variacion genetica fue mas pronunciada en las hembras que en los machos. 
Como en otras aves, esto probablemente se debe a que las hembras son menos filopatricas que los 
machos. El numero de migrantes entre Hungria y Suiza es cercano a un individuo por generation, y 
los valores de Nm calculados por separado para machos y hembras son de 0.525 y aproximadamente 1, 
respectivamente. Es probable que la diferencia en el numero de migrantes se deba nuevamente a la 
mayor filopatria de los machos. La diferenciacion sexual es mas pronunciada en la poblacion suiza que 
en la hungara, y estas poblaciones estan sustancialmente subestructuradas geneticamente. Posibles ra- 
zones para explicar la considerable subestructuracion podrfan ser el comportamiento no migratorio y 
el apareamiento monogamo de la especie, ademas de la separation de las poblaciones estudiadas por 
barreras geograficas, como los Alpes. 

[Traduction del equipo editorial] 


In a number of birds, color polymorphism has 
evolved in allopatry. Under this scenario, a geo¬ 
graphical barrier has separated a population in 
two, a process that facilitates genetic differentia¬ 
tion as gene flow is physically reduced. In some 
cases, during the time of allopatric separation, sto¬ 
chastic processes or natural selection promoted 
the evolution of alternative color morphs for which 
the expression is under genetic control and not 
sensitive to the environment. After the geographi¬ 
cal barrier was removed, the two subpopulations 
could mix in a zone of secondary contact. If hy¬ 
brids are viable, color polymorphism allows re¬ 
searchers to trace back the origin of individuals. 
For example, in the Barn Owl (Tyto alba ) the white 
subspecies, T. a. alba, is found mainly in southern 
Europe, whereas the reddish-brown subspecies, T. 
a. guttata, occupies northern and eastern Europe. 
In central Europe, the two subspecies seemingly 
pair randomly with respect to coloration (Baudvin 
1975, Roulin 1999, Matics et al. 2002), which im¬ 
plies that hybrids may not be selected against. 

Following Voous (1950), these two subspecies 
may have evolved because during the last ice age, 
Barn Owls subsisted in two refugia located in 
southwestern Europe ( T. a. alba) and southeastern 
Europe (T. a. guttata). After the ice retreated, the 
two subspecies invaded Europe via Spain and 
France (T. a. alba) and the Balkans (T. a. guttata) 
to meet in a zone of secondary contact in central 
Europe. The existence of a color polymorphism is 
consistent with the hypothesis that Barn Owl pop¬ 
ulations located in central Europe involve two sub¬ 
species. From a genetic point of view, we expect 
that genetic variation should be more pronounced 
in the zone of hybridization in central Europe than 


to the east or west. The various extant populations 
show differential degrees of introgression suggest¬ 
ed by the phenotypic distribution of individuals 
(Matics et al. 2002, Roulin et al. 2001). In this case, 
one would expect disequilibria of neutral markers 
but this was not tested with our methodology. The 
transition zone appears to be very wide in compar¬ 
ison with other bird species, ranging from west- 
central France to eastern Hungary and northeast¬ 
ern Poland. On the eastern side of the zone, the 
reddish-brown form is more frequent comprising 
84-92% of the individuals in Hungary (Matics and 
Hoffmann 2002) and 90% in southeastern Ger¬ 
many (Schonfeld 1974), whereas the white subspe¬ 
cies prevails on the western side of the zone (Glutz 
von Blotzheim and Bauer 1980; e.g., 75% in cen¬ 
tral France [Baudvin 1975] and 50% in Switzerland 
[Roulin et al. 2001]). The closer the distribution 
of the subspecies is to the 1:1 ratio in a population 
the greater variability of the population due to the 
differences in genetic constitution of the subspe¬ 
cies (Harrison 1993, Arnold 1997). In addition, be¬ 
cause the species is socially monogamous and non- 
migratory, we predict a relatively high genetic 
substructuring of the populations. Because in this 
species, as in other birds (Greenwood 1980), males 
are more philopatric than females, we expect that 
in both Swiss and Hungarian populations females 
are genetically more diverse than males as shown 
in the Black-billed Magpies (Pica hudsonia ; Wang 
and Trost 2001). 

We tested these two predictions by quantifying 
genetic variation using random amplified polymor¬ 
phic DNA (RAPD) technique. This technique is 
based on the amplification of unknown DNA se¬ 
quences using single, short, random oligonucleo- 
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Figure 1. Location of the Barn Owl populations examined for genetic variability. The BIOTA© Association (Pecs, 
Hungary) provided the base map for this figure. 


tide primers. RAPDs provide an unbiased sample 
of DNA variation along the entire genome (Hwang 
et al. 2001) including non-nuclear genomes. The 
sensitivity of RAPDs to population divergence may 
be derived from rapid evolution of non-coding, re¬ 
petitive DNA sequences detected by RAPDs (Plom- 
ion et al. 1995). Therefore, the RAPD technique is 
able to detect variation within and among wild 
populations (Haig et al. 1994, Horn et al. 1996) 
and among sexes within a population (Wang and 
Trost 2001). In birds, mainly species with insular 
distribution (Zwartjes 1999) or living in fragment¬ 
ed habitats (Bouzat 2001) have been analyzed. 
However, this technique has also been used for 
other purposes, including sex determination (Park 
et al. 1997), analysis of wild versus captively-reared 
populations (Bagliacca et al. 1997), and detection 
and eradication of hybrids (Negro et al. 2001). In¬ 
formation on the amount of genetic variation with¬ 
in a species and its distribution within and between 


populations aids conservation planning as well 
(Hwang et al. 2001). 

Study Area 

The study areas were located in Hungary (between 
47°02'N, 17°33'E and 45°50'N, 18°29'E) and Switzerland 
(between 46°56'N, 7°03'E and 46°44'N, 6°41'E; Fig. 1). 
The Swiss study area is located in the middle of the Eu¬ 
ropean distribution of Barn Owls and has an elevation 
between 400-650 m above sea level (masl); the Hungar¬ 
ian study area is near the distribution limit of the species 
with an elevation between 50-500 masl. 

Methods 

Sample Collections. Blood samples from breeding 
Barn Owls were collected between May and August 1998. 
First, the skin was cleaned with ethanol over the point 
where the brachial vein crossed the elbow and insulin 
needles (diameter = 0.4 mm) were used to collect 100- 
|aL blood. Samples were stored on ice in sterile 1.5 mL 
Eppendorf tubes until they were carried to the laboratory 
and frozen at —20°C. The birds were held for at least 5 
min after blood collection and then were released. 
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Table 1. Values of Nei-Li genetic distance among individuals within the groups given. Mean, SE, and Confidence 
Intervals (CT) are estimated by jackknife procedure (H = Hungary, CH = Switzerland). 



Mean ± SE 

95% Cl 

99% Cl 

N 

H females 

0.358 ± 0.006 

0.346-0.370 

0.342-0.374 

13 

H males 

0.261 ± 0.010 

0.241-0.280 

0.235-0.287 

7 

CH females 

0.501 ± 0,021 

0.460-0.541 

0.447-0.554 

17 

CH males 

0.368 ± 0.010 

0.349-0.387 

0.342-0.394 

15 

H pooled 

0.338 ± 0.004 

0.331-0.344 

0.329-0.347 

20 

CH pooled 

0.471 ± 0.008 

0.455-0.487 

0.451-0.492 

32 


DNA Extraction. Standard phenol-chlorophorm-isoa- 
myl alcohol methods (Sambrook et al. 1989) were used 
to isolate total DNA from blood samples. Fifty |xL-blood 
samples were suspended in 200-|xL PBS buffer and cells 
were sedimented. The resulting pellet was suspended in 
extraction buffer containing 20 pg/mL RNAse and was 
incubated at 37°C for 1 hr. Proteinase K treatment was 
then applied (to a final concentration of 100 pg/mL) 
and incubated at 50°C for 3 hr. Samples were extracted 
three times with equal volumes of phenol followed by 
ethanol precipitation. DNA was washed with 70% ethanol 
and resuspended in 100-200 pL sterile-distilled water. 
The concentration of die extracted samples was quanti¬ 
fied using the photometric device GeneQuant (Pharma¬ 
cia, Cambridge, U.K.). 

RAPD PCR Procedure. To find the optimal concentra¬ 
tion, where the reactions gave the most consistent and 
clearest products, DNA was diluted in the range of 10- 
80 ng/pL (by doubling dilutions). In standard experi¬ 
ments 20-ng DNA in 25 pL reaction volume was used in 
die presence of 1.5 mM MgCl 2 . Amplifications were car¬ 
ried out using PTC-150 Minicycler (MJ Research, San 
Francisco, GA U.S.A.) with the following program: after 
a first cycle (2.5 min at 94°C, 1 min at 35°C, and 2 min 
at 72°C), additional 35 cycles were done (40 sec at 94°C, 
40 sec at 35°C, and 1 min at 72°C). The products were 
run on 2% agarose gels. Gels were photographed using 
the digital-gel documentation system BioDocIt (UVP, 
Cambridge, U.K.). Twenty-one different primers were 
tested (QUIAGEN, Budapest, Hungary) and those giving 
the most variable patterns were used for further analyses 
(OPW-17, OPW-08, OPT-20, OPN-05, OPO-05, OPH-14, 
and OPJ-12). 

Scoring and Data Analysis. Fragments were visualized 
by adding 0.1 pg/mL ethidium bromide to the agarose 
gels. Gel photographs were scored for the presence or 
absence of RAPD bands. A pair-wise matrix of the genetic 
distances between individuals was obtained using a eu¬ 
clidean distance measure (Huff et al. 1993), calculated 
from presence or absence data using RAPDistance 
(Armstrong et al. 1994). Components of variance parti¬ 
tioned into within- and between-populations were esti¬ 
mated from this matrix using AMOVA program version 
1.55 (Analysis of Molecular Variance; Excoffier et al. 
1992). The number of permutations for significance test¬ 
ing was set at 1000. Where O sx values differed signifi¬ 
cantly from zero, the number of migrants per generation 
was calculated using Nm = 0.25 (d’sT 1- !) (Wright 1951). 
AMOVA variance components were used as estimates of 


the genetic diversity within and between populations. Be¬ 
cause the genetic distance values were not independent 
data points, we calculated Nei-Li (Nei and Li 1979) ge¬ 
netic distances as well, followed by jackknife analysis 
(Shao and Tu 1995), to find differences between popu¬ 
lations and genders. 

The RAPD technique has some limitations. The most 
important might be that the banding pattern produced 
represents nuclear (both sex-linked and autosomal) and 
mitochondrial genomes. The fact that mitochondrial 
genes are exclusively maternally inherited, whereas au¬ 
tosomal ones are biparental markers, may obscure inter¬ 
pretation of results, particularly in connection with sex¬ 
ual differentiation and consequences of sex-biased 
dispersal (Goudet et al. 2002, Prugnolle and de Meeus 
2002, Vitalis 2002). 

Results 

The amplifications with seven primers produced 
106 reproducible bands in both populations (£ = 
15.1 bands per primer). In the Hungarian popu¬ 
lation, 80 of these bands were present and five of 
them were invariant, whereas in the Swiss samples, 
104 were present and no invariant band was found. 
There were six singletons (bands with one inci¬ 
dent; i.e., occurring only in one individual). 

The 99% confidence intervals of the two sexes 
did not overlap. In both populations the jackknife 
estimates of mean Nei-Li genetic distances of fe¬ 
males were greater than those of males (Table 1). 
The Swiss population showed higher jackknife es¬ 
timate of mean genetic distance than the Hungar¬ 
ian population because there was no overlap in the 
99% confidence intervals (Table 1). 

Because the probability that random distance 
(<fi ST ) was greater than the observed distance was 
less than 0.05 in all cases, the d^g^-values were sig¬ 
nificantly different from zero. The AMOVA indi¬ 
cates 7% of the variance is partitioned between the 
genders within the pooled sample. The Hungarian 
sample gave similar results (7.2%), whereas the 
Swiss population was of higher value (14%; Table 
2). The among-countries differentiation reaches 
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Table 2. Results of the analysis of molecular variance with a between-gender and among-countries arrangement. 



df 

SS 

MS 

Variance 

Component 

(pST 

i 53 

Hungary 







Between genders 

1 

21.15 

21.15 

00.966 (7.2%) 

0.07 

<0.001 

Among individuals/within genders 

18 

222.55 

12.36 

12.364 (92.8%) 



Switzerland 







Between genders 

1 

57.91 

57.91 

2.637 (14%) 

0.14 

<0.001 

Among individuals/within genders 

30 

476.28 

15.88 

15.876 (86%) 



Pooled (all females vs. all males) 







Between gender 

1 

49.30 

49.30 

1.278 (7%) 

0.07 

<0.001 

Among individuals/within gender 

50 

843.50 

16.87 

6.870 (93%) 



Pooled (all Hungarian vs. all Swiss) 







Among countries 

1 

114.92 

114.92 

4.04 (21%) 

0.21 

<0.001 

Among individuals/within countries 

50 

777.89 

15.56 

15.56 (79%) 




a Nonparametric randomization test with 1000 permutations. 


21% (Table 2) and the estimated number of mi¬ 
grant individuals between the Swiss and Hungarian 
populations is 0.96 per generation. The pair-wise 
O sx values gave the same results as the between- 
sexes arrangement (Table 3). Moreover, the 
among-males (0.32) and the among-females values 
(0.20) suggest a higher genetic difference of males 
between populations. From these data we estimat¬ 
ed that 0.53 males and 1.01 females migrate be¬ 
tween populations each generation (Nm values). 

Discussion 

Genetic differentiation among breeding females 
was greater than among breeding males in both 
Hungary and Switzerland. The differential dispers¬ 
al distances of the sexes, males being more philo- 
patric than females, may explain this differentia¬ 
tion. Therefore, males are genetically more similar 

Table 3. Pairwise genetic distances of the populations 
(4> sx between pairs of populations). Above diagonal is the 
probability that random distance (4> sx ) > observed dis¬ 
tance; (number of iterations = 1000; H = Hungary, CH 
= Switzerland). 



CH 

Male 

CH 

Female 

H 

Mail 

H 

Female 

CH male 

— 

0 

0 

0 

CH female 

0.14 

— 

0 

0 

H male 

0.32 

0.29 

— 

0 

H female 

0.25 

0.20 

0.07 

— 


to each other within a population. This scenario, 
based on sex-specific dispersal, may also explain 
the greater among-males differentiation between 
populations (about twice as many females as males 
appear to emigrate). Because males do not emi¬ 
grate as much as females, they preserve the genetic 
features of their population to a greater degree. 
This is further supported by the fact that the esti¬ 
mate of migrants per generation (based on genetic 
data) gave the same result as the ringing data: 
greater male philopatry (Taylor 1994). Previous 
studies suggested that, from Hungary, a greater 
percentage of individuals move toward central Eu¬ 
rope than in the opposite direction (Matics 2003); 
i.e., the Hungarian population was more of a 
“source” than a “sink” population. As a conse¬ 
quence, gender- and population-differentiation is 
greater in Switzerland than in Hungary, because 
the exchange of individuals is guttata- and female- 
hiased. These results were concordant with the fact 
that both phenotypic and genotypic variability of 
individuals were greater in the middle of a transi¬ 
tion zone than on the edges (Arnold 1997, Roulin 
2003). This proposal could be tested with data 
from other localities such as from western Euro¬ 
pean Barn Owl populations. 

The between-gender differentiation of the spe¬ 
cies detected using RAPDs (7%, 14%) seem to be 
disproportionally high. Using the random priming 
technique the sexual differentiation was expected 
to be between 1-2%, as the Barn Owl has 92 chro- 
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mosomes (Belterman and De Boer 1984). This 
higher value of gender differentiation could be 
due to the relatively large size of sexual chromo¬ 
somes in birds (Stevens 1991). The sexual differ¬ 
entiation of autosomal markers caused by sex- 
biased dispersal is instantaneous because the next 
generation receives a random set of alleles from 
both parents. On the other hand, sex-chromosom¬ 
al markers preserve the differentiation for longer 
time and when gene flow occurs continuously, this 
differentiation could be detected permanently. Dif¬ 
ferentiation detected in this study therefore may 
be associated with markers located on sex-chro¬ 
mosomes. 

Although only two populations were analyzed, 
the O st value of 0.21 indicated a substantial genet¬ 
ic substructuring among our study populations. 
The results of another RAPD study gave <P ST values 
of 0.048 and 0.103 for island species (two popula¬ 
tions of both the Puerto Rican Vireo [ Vireo latimeri ] 
and Jamaican Vireo [V. modestus ], respectively), 
and 0.015 for a migratory continental species 
(three populations of the White-eyed Vireo [ V. gri- 
seus]; Zwartjes 1999). For the Greater Rhea (Rhea 
americana) a 4> sx value of 0.0637 was found among 
four wild and a captive population (Bouzat 2001), 
which is a low value given that this species is flight¬ 
less. The special mating system of the Greater Rhea 
could play an important role in producing this low 
genetic substructuring. The male rhea establishes 
a territory and builds a nest. He will then attract a 
group of about three to six females with whom he 
mates and they lay ca. 20-30 eggs in his nest. While 
the females go off to mate with other males, the 
male will incubate the eggs and look after the 
chicks on his own. The relatively strong substruc¬ 
turing of the Barn Owl populations could be ex¬ 
plained by at least three factors: (1) nonmigratory 
behaviour, (2) the socially monogamous mating 
system of the species, and (3) the presence of a 
geographic barrier (Alps) between the populations 
analyzed. The Nm value of roughly one is the min¬ 
imum amount of gene flow that prevents differ¬ 
entiation at neutral loci among populations by ge¬ 
netic drift (Wright 1931). 

The conversion of F ST -related values into Nm is 
problematic for several reasons, including that it is 
based on isolation-by-distance models. In this 
study, we could not test for correlations between 
genetic and geographic distance. Furthermore, 
population size and dispersal rates are not constant 
over time and space and assumptions of demo¬ 


graphic and genetic equilibrium and uniformity 
are unrealistic (Whitlock and McCauley 1999). 
Many assumptions of the models used are violated, 
so that results cannot be interpreted as direct mea¬ 
sures. In addition, other evolutionary forces con¬ 
tribute in establishing differentiation (Bossart and 
Prowell 1998). Finally, we suggest that drift should 
play an important role in the microevolution of the 
Barn Owl as this species of tropical origin (Voous 
1988) goes frequently through bottlenecks in the 
suboptimal European area in hard winters (Taylor 
1994). 
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